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ABSTRACT 


Paiperramewer=plane method, baséd upon a quasi-frozen 
assumption, for the design of electrical filters containing 
Slowly varying nonlinear elements is presented together 
with examples»sof its application vo the lowpass, highpass, 
bandpass and band-reject filter types. The nonlinear mag- 
nitude and frequency characteristics are presented in the 
form of parameter-plane graphs which allow value judgements 
to be made concerning the choice of nonlinear elements used 
in the filter design. The computer programs used to gener- 


ate the graphs are also included. 
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I. INTRODUCTION 


A filter is a frequency-Sensiftive unit which discrimin- 
ates between desired and undesired frequencies. Ideally 
the filter would have a unity-gain passband and a zero- 
gain stopband with no transition region between the bands. 
Also, its phase shift would be linear within the passband. 
such filters do not exist today. The ideal fi ieee ele 
optimum case and serves as a Standard of compareegi 107 
all other filters. Filter design usually stresses either 
the magnitude response or the phase response depending on 
the intended use of the filter. 

Filters may be grouped into four distinct types based 
on their frequency discriminatiom: lowpass) Weems 
bandpass, and band-reject. While classical design used 
passive, lumped, linear elements and image-parameter methods 
to realize these filter types, modern design uses the poly- 
nomial method of network synthesis and considers all types 
of elements in both their active and passive forms. The 
addition of nonlinear elements in these filters allows 
a closer approach to the ideal response than possible if 
the elements were linear. At the same time however, the 
design problem becomes more complicated. 

A parameter-plane method for the design of filter 
circuits is presented in this thesis. The method provides 


a relatively simple, easy-to-use way of incorporating 


x 






nonlinear elements into the design of electrical filter 
——-> > =p ‘ 


CiIeeuLES. 
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IL. THE PARAMETER PLANE MEAOp 


The mathematical analysis of a linear dynamic system 
usually evolves into the evaluation of the roots of a poly- 
nomial formed from tne differential equation descr vm om 
the system. The root locations of this characteristic 
equation are then interpreted with respect to the system 
response cNaracteristics. Any Change in “ally eaomee 
parameter value, however, changes the root locations. S@me 
changes will improve the circuit response while others will 
degrade it. 

Mitrovic, Thaler, Siljak and others developed analysis 
methods which consider such parameter value variations [1]. 
The original method of Mitrovic was rather limited since 
only the Last two coecificients of the characteruerune 
equation were allowed to vary. Lines of constant damping 
ratio (t) were mapped from the s-plane into the plane of 
Che variable coefficients using the characteristic equation 
as the mapping function. These transformed curves were 
Chen used to determine the variable coefficient values 
required to place the characteristic €qUamiien Voetserar 
desired locations. 

Generalization and extension of Mitroviec Ss coeriicieng: 
plane by Siljak and Hollister [2] led to the concept of the 
a,8 parameter plane. This method considers the variable 
coefficients of the characteristic equation as nonlinear 


functions of two system parameters, oa and 8. 
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F(s) ,s = O, a. = f.(a,B). 


v 


ll 
wm bd 
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Curve families of constant_damping raviio pee undamped nat- 
ural frequency are commonly mapped into the a,8 parameter 
plane and used to infer system characteristics. 

In the case of a system transfer function, however, the 
Chareacuericule equation is formed from the -denomna por 
tic aunertatoreis mot includedsin the»formmul atdion. Jipeve= 
fore, full frequency-response information cannot be gained 
from using the characteristic equation in the parameter 
plane method. Thaler and Thompson developed a form of the 
transfer function which could be used in the parameter-plane 
method [3]. This transfer function form has been used by 
Glavis [4] and Staples [5] to make frequency-response plots 
and analyses of systems in which the coefficients contained 
iba? Cali angie Ons. Of 0. 81 sb. 

The derivation of the transfer function form developed 
in Ref. 3 follows: 

Let the system transfer function be given as 


N(s) 


ms) = Dts) 


where NS and D)(s) “are each of the’ form 


Nts) = as. 


i 
- 


dy? 


and the terms are nonlinear function of two parameters, 
a and 6p. The squared magnitude function can be obtained 


from 


2 . ; 
nm =dhigjiw) Dhegeve 
If the polynomials N(s) and D(s) are separated into their 


even and odd parts, 


N.(s) + N_(s) 
Gee O 
MSs) = 5°(s) + DCs) ° 
MW is) ee) 
=, 0 
is D.(s) ss Dos) 


N (s)° SON (s)- 
and T(s)-T(-s) = 2 _, 
D.(s)* = D(s) 


From this it may be seen that both the numerator and denom- 


inator are even polynomials. When s = jw both numerator 
and denominator are even in w. Thus 
si ; 
25 A,wo* 
me 7 * 
a Bw 
k=0 
where A. = f,(1,0,8,08,0°,8°,0°B,08°,...) 


and similarly for the B 


ee EUS 


There are now four parameters present, alpha, beta, 


magnitude and omega. Any one may be expressed in terms of 


1s: 


~*~ 


the others and curves drawn depicting the various relation- 
ships among the parameters as defined by the transfer 
function. Phase versus omega may also be examined by 
rationalizing T(jw) and using the arctan relationship to 
solve fer the angie. 

Although the original motivation for the parameter plane 
was to analyze linear dynamic systems, the method may also 
be used to analyze nonlinear systems. Consider a linear 
system. One of its characteristics is that a sinusoidal 
input signal of a given frequency will be repeated through- 
out the system as a sinusoidal signal of the same frequency. 
LitcmomiCwmaiateaws Crue in a clreuit coOmpaining nonlinear 
elements. In most cases though, the internal signals are 
at least periodic and have the same period as the input 
signal. These signals may therefore be represented by a 
Fourier series. Thus a frequency-response transfer function 
may be defined relating the fundamental component of the 
output signal to the input signal and may be used in the 
parameter-plane analysis of a nonlinear system. The non- 
linearity is considered a parameter whose characteristic is 
Gefined by a Gens of POoINntS In a padwameger plane e.g&. ,; 
alpha versus transfer-function magnitude. The function 
a(v) is assumed to vary slowly enough that it may be con- 
sidered Pitenn sen! or fixed over the frequency range of 
nee re ann 

Purtvher, this parameter plane analysis logic may be 


reversed and the method used as a tool in the design of a 
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system. Thus a desired system response may be specified 
and the «2,8 parameter variations necessary to produce this 
response determined by the method. 

The parameter-plane method would be far too laborious 
if done by hand. Therefore, FORTRAN IVecomputer pregram 
have been developed by the author which allow the frequency- 
response form of the parameter plane method to be used in 
the design of filters with nonlinear elements. The program 


details are contained in Appendix A. 
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Filters with nonlinear elements may be designed by 
reference to parameter-plane plots of the voltage transfer 
Purceror Of ele 11iter. The design may ofpneeive the mag- 
nitude response, the phase response or both. The nonlinear 
elements are represented as parameters alpha and/or beta 
and used in conjunction with the parameters of magnitude 
and of normalized radian frequency, omega. A computer 
generates selected plots of any one parameter versus another 
with a third parameter held constant and the fourth incre- 
mented to form a family of curves. A fifth parameter, 
phase, may be plotted versus omega for family parameters 
alpha and beta. The desired filter response is noted on 
these plots and the required characteristics of a particular 
nonlinear element determined. 

reo CclveMmespecifications @ particular type of filter 
Civico senmeanGdsthe transier functlon determined: 

For simplicity let the nonlinear element be represented by 
the alpha parameter and be a function of output magnitude. 
Also assume that in this case the magnitude response is 
important while the phase response is not. 

The computer is told to make a Bode plot of the linear 
response of the filter for each of a family of alpha values 
and the response required by the specifications isS over- 


plotted by hand onto the same graph (see Figure 1). Each 
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alpha-omega coordinate point along the specification curve 
indicates what the value of alpha must be in order to 
satisfy the specifications. These coordinate points are 
then picked off and transferred by hand to a computer- 
generated parameter-plane graph of alpha versus magnitude 
where beta has been held constant and omega incremented to 
form the curve family (see Figure 2). The locus of hand- 
plotted points defines the required charaecverrscvicenol 

the particular nonlinear element with respect to the circuit 
output magnitude. If the nonlinear element were a function 
of frequency, the coordinate points would have been hand 
transferred to a computer-generated parameter-plane graph 

of alpha versus omega with magnitude as the family parameter. 
This would determine the necessary frequency characteristics 
of the nonlinear element. 

Filters with two or more nonlinear elements may require 
plotting points in the alpha-beta parameter plane (magnitude 
held constant and omega incremented or vice versa) to 
determine possible value trade-offs between the nonlinear 
elements. The curve plotted on the graph of parameter ver- 
sus magnitude or omega would still be used to define the 
nonlinear element characteristics. These same procedures 
could be followed with the response emphasis on phase rather 
than on magnitude. 

fin additional consideration is necessary." =the trageier 
function is the ratio of one voltage to another, E2/El. 


Consequently, the plotted response is that of E2 with El 
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set to onemwvolt. Thus in realizing the filter Gece 
nonlinear element must be made a function of E2/El and not 
a function of E2 alone. <A circuit to _ perfoirnpieiis Soieee een 
may be omitted if the nonlinear element is a function of 
frequency or when simulation sensitivity checks indicate 
that the filter response will be within acceptable limits 
for the predicted range of input values. 

Although the filters considered here are shown in a 
lumped R - L - C configuration, the R, L and C should be 
looked upon as representing voltage-current relationships 
rather than actual physical components. The lumped R-L-C 
could be transformed into distributed values and/or L - C 
interchanges made through the use of active networks such 
as gyrators, mutators, rotators, ete. Time 1nvami aieeee 
assumed however. 

The steps in the design of a filter are: 

1. Form the transfer function, T(s), of the desired 
filter type designating the nonlinear elements as alpha, 
beta. 

2. Determine which computer-generated plots will be 
meeded: 

3. Encode data and use the computer programs described 
in the appendix to generate the plots. 

4, Select the Bode and/or phase plot as appropriate 


and draw in by hand the specified response curve. 


. 


5. Hand transfer the coordinate points along the res- 
ponse curve to the computer-generated parameter-versus- 
magnitude Or parametver-versus—omega prapns. 

6. Determine the characteristics of the nonlinear 
element trom the loeus of points plotted in step five above. 

(oaeiechweene Cdesien Dy Computer Simulation using 
several voltage input values to determine the effect of a 
Chawieiime INpul. 

8. If the nonlinear element is a function of magnitude, 
decide whether or not the divider circuit is necessary. 

9. Frequency shift the component values from the 
normalized radian frequency to the frequency of interest. 
Impedance shifting may also be used. 

In the case of two or more nonlinear elements, step five 
will include the alpha-versus-beta plot for possible | 


parameter value trade-offs. 


A. FILTERS WITH ONE NONLINEAR ELEMENT 

The following filter examples cover the four frequency 
types, lowpass, highpass, bandpass and band-reject. Each 
of these basic filters is made up of three elements, one R, 
one L and one C with the nonlinear element being a function 
of the output magnitude. The R element has been chosen as 
the nonlinear element since it may be readily associated 
with a field-effect transistor. The L or the C could just 


as well have been chosen. 
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ee Lowpass Filter 


Figure 3 illtistrates Gne possible lempass eer 
conTrgunratbton. Since frequency scaling willbe employed Bs 
the last design step, the L and C values are normalized to 
one. The parameter alpha replaces the nonlinear clement Re. 


Thus Ghe circuit transfer functicmiis 


i 
lore + Ss 


T(s) = 5 
The necessary. computer plots, for the desiegm@ Gf jigs 
circuit are the magnitude-versus-omega (Bode) and alpha- 


wersuSs-magpnitude plots. 





Figure 3. Lowpass Filter 


Figure 4 is the Bode plot for the circuit. The alpha 
parameter curves on the plot indicate the output magnitude 
of the filter for the normalized radian frequency range of 
O.1 to 10.0 radians when alpha is fixed-at a parereular 
parameter value. Thus if alpha were set equal to 1.0, the 
output-versus-frequency response of the filter would be 
that indicated by the curve for alpha equal to 1.0. There- 


fore, this graph illustrates the response of aW@inear tite 


on 


mee | 





WH 0.90 


FIGURE 4. 
Ro»Ae-NICHGLS, RODE PLOT. PARANETER-ALPHS 
LGHPASS FILTER. L-C=1.0. R-=ALPHA 


Ze 


for several values of one of its elementsm Raper 
specifications are indicated by the hee Ixie dragon 
the graph by hand. In this case it is assumed to be the 
magnitude response of an ideal filter. )UThewzeiro cain 
the stopband of the ideal filter is represented by a@ Time 
at the minus 60.0-dB level. 

The coordinate points of the intersection of the alpha 
parameter curves and the desired response line in Figure 4 
are picked off and transferred by hand to the alpha—-versus— 
magnitude parameter=plane plot in step fiwe of the @esmen 
procedure. This plot is Smewn ine@igure 5. Sieve 
line is the locus of points transferred. It indicates the 
desired characteristics of the nonlinear element when The 
Value of the element is a funetaonsofethestiaiterwsouinan: 
magnitude. 

Mre initial Wwalwe of alpha for This plot 3s lesser, 
orveqcual tome .@rat “an omega of 0.1. The required value 
of alpha at tcutoffr is 1.0.% For a flat, zero=@b ros ose 
between these two omega points, required alpha values are 
1.5 at ah omega of O.2 and 1.35 at an omega equal to By. 
The plot also indicates that an alpha value of 1.0 would 
deviate from the ideal response by +1.0 dB at omega equal 
to 0.5eand less thaneabout +1.2 dB betweem 0.5 anm@ecutorr. 
If this is tolerable, alpha may be fixed at 1.0 for values 
of omega less than 1.0 and the design problem simpamefied. 

In the region where omega is equal to 1.0 the alpha- 


mMapnitude relationship is logarithmic. Themreivereeurpur 


ZS 
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FIGURE 5. 
ReA.NICHOLSsALPHA US MAGNITUDES DB} .PARRN=ONEGAR 
LGWPASS FILTER» L=C=1.0. R=ALPHA 
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0.00 
MAGE DB} 


is 0.0 dB at alpha equal to 1.0, -60.0 dB at alpha equal 
to 1000 and -80.0 dB at an aSlphasof 10000. “Ties 
60.0-dB level requires alpha to have a three-decade range 
in value. The minus 80.0-dB level is attainable with the 
four-decade range of Pilea toes transistors. Alpha may 
remain at this last value for frequencies abovewiaGe Tt 
this isso, the output leveli will falloff at ae2zeeae 
decade rate. 

At this point consideration should be Given tome 
realization of the circuit.@ Although theoretregi sa) 
function may be generated, the economics of the situation 
usually limits the degree of circuit complexity allowed 
and causes simplifying approximations to be investiigaren. 
The fixing of alpha at 1.0 for the frequency region belem 
1.0 and to a fixed value for the region above 1.0 would be 
Simplefileations in Ghesease of the lowpasteiilver. “Tire 
drop in output magnitude below 0.0 dB at cutoff when alpha 
Ws equal to 1.0 may be used@to initiate the reqmired 
parameter value change. The change in K could come fiom 
the filter transfer=funetion output voltage Deine ewer ed 
from a fixed reference level and the difference applied To 
the gate of a field-effect transistor. 

Figure 6 is the phase-plot companion to the Bode plot 
of Figure 4. The dark line indicates the expected phase 
response of the filter when alpha equals 1.0 at values of 


omega below 1.0 and equals 10,000 at omega values above 1.0. 
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FIGURE 6. 
Re R»NICHOLS, PHASE PLOT. PARANETER=RLPHA 
LOWPASS FILTER. L=C=1.0. R=ALPHA 
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The. phase. fall-offiin. the latter region is very slow. 
Designing for the ideal magnitude response does not yield 
a linear phase filter. 

If alpha were a functi@n of frequency, the wou 
alpha versus omega, Figure 7, would replace the alpha- 
versus-magnitude plot. Again, the heavy line is the locus 
of points transferred from the Bode plot of Figure 4. The 
O.0-dB curve defines the exact values of alpha necessary 
for a completely flat magnitude response in the passband 
region. The area between the 0.0 and +1.0 @B curves indi- 
cates the latitude in alpha values available for frequencies 
below 0.2 rad#ans. 

The c@mputer simulation and input SensiUiViiy gene, 
results are depicted in Figure 8 for the output magnitude 
and in Figure 9 for the phase. A 1.0-dB region of uncer- 
tainty or lag is assumed fer the output level av wWikstem Uie 
parameter value change is initiated. Alpha is fixed at 
1.0 for output values above the minus 1.0-dB level and 
allowed to go as far as 10,000 whenever the output drops 
below -1.0 dB. The curves répresent the circuit» response 
when the divider circuit is omitted and the input voltage 
Varies ar@und 1.0. 

co? H@ehpass Filter 

Fi@eure 10 illustrates a possible@mienpess filve. 
configuration. As in the Mowpass filter circuivy 1 Geeme 
are normalized to one and alpha represents the nonlinear R 


in the» transfer function 
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5° 
TCs) = —_ 3 
L@reres + s 


The necessary @domputer plots for the design Bf this jieer 


are the Bode and the alpha-versus-magnitude plots 


ty 


Figeure 104 @ Htehvasse Filter 


Figure ll is the Bode @plot for @his cirGhit. Sitecewee 
circult configuration is related to that of the lowpass 
filter by an interchange of frequency-sensitive elements, 
the Bode plot isfla reflection of, the Wowsess Pilver Girevit 
about the omega equal to 1.0 axis. This relationship is 
further brought out by the alpha-versus-magnitude plot of 
Figure l2. It is exactly the same as Figure 5 except for 
reversed omega values. 

Cireuilt realgzavion is agdso related tO UiawOol a. 
pass filter. InjJthis case the parameter alpha coudtd@be 
fixed at 10,000 in the omega region below 1.0 and allowed 
to decrease to 140 in the region above 1.0. The filter 
transfer-function output would be added to a reference level 
and the sum applied to the gate of a field-effect transistor. 
The change in alpha would be initiated when the output mag- 


nitude exceeded a given value, say -80.0 dB. 
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3. Bandpass Filter 
The three-element circuit configuration for a band- 
pass filter is shown in Figure 13. Again L and C are set 


to one and R replaced by alpha. The transfer function is 


T(s) gee 
lL + @s + 5 


The Bode plot and the alpha-versus-magnitude plot are the 


computer plots used in the filter design. 


LOA 





Figere 13. ¥Bbandpass Pailter 


The heavy line drawn on the Bode plot of Figure 14 
specifies a 1.0-omega-wide bandpass symmetrical about the 
pelnt whergememesea ecGuals 1.C0@and a stopdand level of 
-60.0 GB. 

Pigures 15 shows the specifications in terms of the 
nonlinear element versus magnitude. The desired alpha 
characteristics in this case are more complex than in the 
lowpass case. A small change in alpha causes a large 
output change when alpha, itself, is small. When alpha is 
large, which is the case in the passband, large parmeter 


value changes are necessary to affect the output. This 


34 





S809 0: a 
MAGE DB) 
FIGURE 14. 
R.A-NICHGLS, BODE PLOT. PARAMETER=RAWPHR 


BANDPASS FILTER» L=C=1.0. R=ALPHA 


20.90 


SD 


Czy 


Won 
04 
08 


V2 





U1 
<0. 00 





Fe eee 
B.A NICHOL SPRL PRHOUS e 
BRHNUFRSS F | re Be Pee AA. R=HL PHH 


HAGONT TUBE! OB), RARRM=QGHEGA 


36 


would indicate a relatively unstable attenuation level in 
the stopband region but a stable level in the passband. 

The width of the passband is limited by the minimum 
attenuation required in the stopband. The wider the pass- 
pand, the higher the maximum alpha value mist be for a Bat 
response. The range of values available with the chosen 
element would then indicate the point where the stopband 
attenuation is too low to meet specificati®ns and thus 
limit the width of the passband. On the other hand, the 
filter passband could be made very narrow. 

Realization of the filter can not be done.on..the.basis 
of magnitude alone. Some sort of frequency or phase 
sensing device is needed to determine the region of opera- 
Tion, either above or below omega equal to 1.0. Given a 
frequency sensor the circuit could be formed by letting the 
magnitude response level crossing initiate the alpha 
parameter value changes to define the end points of the 
Bassband.. Thig..ais the same typesof operati Onsaas@gusedudm 
the lowpass and highpass filters. 

The phase plot is shown in Figure 16. The zero phase- 
shift region indicates the passband where alpha values are 
large. The phase below the passband would be 90 degrees 
leading and above the passband 90 degrees lagging. 

4. Band-Reject Filter 

A three-element band-reject filter is illustrated 
jm Figure 17. The transfer™function when L and © are sserv 


to one and R replaced by alpha is 
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The design of the filter calls for Bode and#alpha@we rame= 


magnitude pleus. 





Figure 17. Band-Reject Filter 


the specifications are indicated on the Bode plot, 
Figure 18, and the alpha-versus-magnitude plot, Figure 
bo. The phase plot is shown in Pigure COFMEhe piers 
confirm the expected, that the band-reject filter is the 
inverse of the bandpass filter. 

As in the case of the bandpass filter, the band- 
reject filter can not be realized by using magnitude res-— 
ponse alone. It too must have a frequency @ phase Seimiar 


to indicate the region of operation. 


B. “FPLTERS WITH MORE» THAN ONE NONLINEAR ELEMENT 
More complex filters may be formed using more than one 


nonlinear element. They may be built up by combining filter 
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sections and/or by using several nonlinear elements in the 
Same SECTiLoOn. 
is Concatemand Tig. Mal Ceres oc iene 
Additional attenuation of undesired frequencies may 
be obtained by concatenating filter sections. The sections 
are designed as individual units and then coupled together 
by bmi.teae amp lim Gees . 

A considerable improvement over the single-section 
design of a bandpass filter may be made by concatenating 
lowpass and highpass sections. The passband end points 
are set in the frequency shifting step by the amount of 
HP-LP cutoff frequency overlapping used. A band-reject 
filter may be designed by inverting the bandpass filter 
SCuUtpDUL and adding it to the input Sigal. 

owe, CAS CeGimewraiter Socrzons 

There’ are two types Of cascaded filvers. Seice wee 
type consists of those in which the nonlinear elememy 
each Section 1s made a function of the transter ragworm 
the total filter.” This type also includes UmoOse Bava. 
Weich the nonlinearity is a function of irrequency gine 
this parameter, or at least its first Fourier approximation, 
me commen Co Che fitter as av-whole. 

Figure 2l shows an example Of "this ifarst "type wor 
cascaded filter. Two lowpass filter Sections areseaseae 
wilvuhn both Rl and R2 functrems ol vwne total Prrver cue 


magnitude. Rl equals k times R2. Setting L, C and k to 


= 


one and replacing the R's with alpha gives a transfer func- 
aor: Jost 


1 


TCs) fib = ela ae a 
2 Ose Ot Sables rs pe + s 





Figure el. Cascaded Filter sections, Type l. 


This type of cascaded filter is designed as a single 
filter since the transfer function contains only one param- 
eter, alpha and it is a function of magnitude. The Bode 
MeLoU mor Chissseircuitois shown in eure 22. the alpha- 
versus-magnitude plot in Figure 23 and the phase plot in 
Figure 24. The Ste 3 are assumed to be the same 
as those for the single-sections lowpass filter,except 
that the stopband attenuation is now set to -120.0 dB. 

In this filter design there is no simple approximation 
which would give a reasonably flat response below cutoff. 
The alpha value must change smoothly from 5.0 at omega equal 
to 0.1 to-0.1 at omega equal to 1.0. The relationship of 
alpha to magnitude in the omega equal to 1.0 cutoff region 


is shown See 25, Tes Is snotea ssi mo lie, 1 oes risiiemec 
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relationship as was the case with the single-section low- 
pass 1 ievex’. 

This type-l cascaded, lowpass filter has an advantage 
over the two-section concatenated filter in that it does 
not need the buffer @mplifier; but, on the other hand,’ it 
does require a more complex parameter characteristic in 
une passdb ang. 

The second type of cascaded filter is illustrated in 
hieure ¢cS.™ Gr consists of@those filters, formed by the 
direct addition of sections. In forming the transfer 
function L and C are set to one, Rl iS replaced by alpha 
and Re by beta. The transfer function for the total filter 
is 


1 


T(s) — H3/E = >... . a, ee 
(at+B)s + (2+0a8)s~ + (atB)s- +s 


Re is a function of E3/E2.% The transfer function for ythis 


second section is 


Tis) = E3/E2 = —__§ Ae 
| +aeosS Bt 2 


E | i prezy hi 39\E3 


Figure 25. Cascaded Filter Sections, Type 2 
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In Gaddisgson, Rl tS,ae function of E2/E1 sotthat 


T(s) = E2/El = 1+ Bs + s° : 
(atB)s + (2+a8)s° + (atB)s> + 5! 

hebpode plot for the E3/El transfer fumetion sie no a 
in Figure 26. The parameter is alpha with beta set to 1.0. 
The@ specification curve is drawn in as in the previous 
examples with the stopband gain level at -120.0 @B. Duewo 
the alpha/beta symmetry in the transfer function, Figure 26 
is also valid for beta as the parameter with alpha set to 
1. @. 

Rl is @ function of an internal voltage-tranmeier rama, 
K2/E1l. The method of design in this case may be illustrated 
by assuming that R2 is a fixed value and that Kis theeanly 
nonlinear element in the circuit. As before, the alpha-omega 
cCoeordinatelpoints of the desigred Pesponse Gurver are: we eed 
off from the Bode plot and transferred to the alpha=versus- 
Magnitude plot illustrated in Figure @/ and ghe Heaumaw 
characteristics determined. The only differemce@nere as 
that the magnitude used is the E2/E1 voltage ratijfo rather 
than the output magnitude of the E3/E1l Transfemencera. 

In this filter, however, R2 is alsow™m@nlinc@ee so Pie 
plot of R2 versus E3/E2 magnitude must be considered. Mis 
plot is the same as that of the single-section lGwpass 
filter in Figure 5. Therefore, realization of the passband 
region may be simplified by setting Re at 1.0 and causing Rl 


to follew the Gharacterishae curve gin Preure 2 / es Biayvalives 
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of alpha to be used in the passbam@@mresion are determined 
from Figure 28. This figure is a plot of alpha versus 
beta for passband frequencies where the filter output 
magnitude is 0.0 dB. 

In the cutoff region Rl and R¢e change from their pass- 
band low values to their stopband high values. R2 
follows the omega equal to 1.0 logarithmic curve to the minus 
80.0-dB level at beta equal to 10,000. A four-decade range 
Rl element goes from 0.1 at 0.0 dB to 1000 at -60.0 dB. In 
so doing it follows the hand-drawn curve in Figure 29. This 
curve is a composite of the line for omega equal to 1.0 in 
Figure 27 where R2 equals 1.0 and that of Figure 29 where 
R2 equals 10,000, The chosen values are determined from 
the E2/E1 alpha-versus-beta plot of Figure 30 where omega 
is fixed at 1.0 and output magnitude is the curve family 
parameter. 

In the stopband region Rl may be left at 1000 and Re 
at 10,000. The resultant attenuation at cutoff is -140.0 
€B from which it falls off at 40 dB/decade. Figure 31 is 
the Bode plot of E3/El1 for the stopband region above omega 
equal tag .0. 

Pigu@mmee 32 illustrates the phase response for passband 
values of alpha and beta while Figure 33 shows the phase 
response for stopband values. 

Figure 34 and 35 show what the output magnitude and 
phase response would have been if two lowpass circuits 


had been added together directly. The response is not the 
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same as would have been obtained by concatenating the two 
sections. The concatenated filter is easier to realize 
than a type-two cascaded filter and has a greater stopband 
eltenwataen. 

3. Two Nonlinear Elements in the Same Section 

A third configuration for two nonlinear elements 
may be obtained by lettime both the K and the © in the@aasic 
lowpass filter be nonlinear elements and functions of the 
filter transfer ratio. The design procedure is the same 
as that of the previous sections. Bode plots are made for 
each of the nonlinear elements with the other set to a best 
guess value. From the specifications and the Bode plot 
information minimum and maximum values for the elements are 
chosen to use in a design trial. The information is 
transferred to parameter-versus-magnitude plots and the 
nonlinear element characteristics determined. The inter- 
dependence of R and C requires a composite curve to be used 
in defining these characteristics. 

Figure 36 is the alpha-versus-beta plot for omega equal 
to 1.0. Alpha represents the nonlinear R which is assumed 
to have a four-decade range of values. Beta represents the 
nonlinear C which is assumed to have a two-decade value 
range. The specifications are those of a lowpass filter. 
The heavy line in Figure 36 indicates the chosen alpha/beta 
relationship in the filter cutoff region. Figures 3/meand 
38 are the alpha- and beta-versus-magnitude plots with the 


element characteristic curves drawn in. Computer simulation 
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resultspifor this design ane shown in Figure 39 foF output 
magnitude and in Figure 40 for the phase. Alpha and beta 
are set to 1.0 when the output magnitude is greater than 
-l.0 GB and allowed to go to 10,000 and 1004.*respectively, 
when the magnitude drops below -1.0 dB. As before, the 
—1.0 dB, level’ is used to représent a lag in changeover fom 
the change@initiation at 0.0 a@B. The,cutoff attenuation 


is 120.0 dB while the phase response is nonlinear. 


G@ DESIGNING) POR PHASE RESPONSE 

When thei filter specifications concern phase Prespon te 
rather than output magnitude the desired response is 
drawn in on the phase plot rather than on the magnitude- 
versus-omega plot. The alpha-omega points are still 
transferred to the alpha-versus-magnitude plot for defini- 
won of the dharacteristics of Ehe nonlinear eclenmne wae 
igs is @& Pimetion of the filter output transt#er fumice iene 
Where the magnitude of the output is of interest the same 
alpha-omega points may also be transferred to the magnitude- 
versus-omega plot for an indication of the magnitude response 
corresponding to the filter phase specifications. 

Figure 41 is the phase plot for the basic lowpass 
falter of Figure. 3... The heavy ine represents 4 Bpe@iiica— 
tion calling for a linear 180-degree phase response between 
omega equal to 0.2 and 1.8. The alpha-versus-magnitude 


plot is shown in Figure 42. The magnitude characteristic 
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FIGURE 41. 
R.A.NICHGLS>»PHASE US GNEGA.PARAN=ALPHA=R 


PHASE DESIGN HITH LOWPASS FILTER. L=C=1.0 
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alpha is fairlyetlal term Paewiewer aie of The pigeon e mioy 
range but falls off in the upper half. This may also@ee 
seen in the magnitude-versus-omega plot of Figure 43. 

The basic lowpass filter may be designed for a 180- 
degree linear phase shift if the output magnitude is not 
part of the specifications. If, however, a flat magnitude 
response is also important, the maximum phase shift would 


be limited to 90 degrees for this frequency range. 
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IV. CONCLUSIONS 


Privers with nonlinear elements may be designed using 
the frequency response form of the parameter plane method. 
The method is a graphical one which might be thought of as 
a counterpart to the describing-function method used in 
analytical design. Parameter-plane plots furnish the mag- 
nitude and frequency characteristics of the nonlinear 
element so that the filter may be designed to meet specifi- 
cations. This method has the advantage of presenting the 
destener with a picture of Whe possdble circuit condi@ions 
and thereby allowing value judgements to be made concerning 
tilesparvclcurar response, chosen. 

On the other hand, this form of filter design has the 
disadvantage common to all graphical methods. As the num- 
ber of variables increases so does the number of graphs 
needed and the rate of growth of numbers of graphs is 
mueh faster than that for the number of vawiables. he 
point is finally reached where the number of graphs needed 
surpasses the ability of the designer to interpret and 
integrate them. At this point a hybrid form of graphical- 
enieusytigge Moe be Used. 

The basic building block for filters with nonlinear 
elements is the lowpass filter. With a nonlinear R having 
a four-decade range of values, it will provide a stopband 
attenuation of 89 dB and a passband with a maximum rise of 


less than 1.5 dB. Adding a nonlinear L or C with a 


(2 


two-decade range of values will increase the stopband 
attenuation to 120 dB. Multiple filter sections may be 
concatenated to reach even higher levels of attenwation. 
If the L and C elements are interchanged, the lowpass 
filter is transformed into a highpass filtem A bandpass 
filter may be formed by concatenating low- and highpass 
sections. The band-reject filter may be genemated fmem the 
bandpass form. The lowpass filter may also be designed 
for a linear phase response. A linear 180-degree phase 
response is available for cases where a flat output maeg— 
nitude is not necessary, or a linear 90-degree response 
where a reasonably flat output magnitude is specified. 
several areas exist for future investigation. The 
first one is that of the physical realization of these non- 
Menear filters, primarily the lowpass filter. Another as 
tae investigation into the response of the [llter ee. 
complex wave shapes. A third area is the extension of the 
mevmoed to include more complicated Circuit Configurar een. 
where the question of the optimum mix of nonlinear elements 
Goes not have an apparent answer. In addition, the use of 
feedback and feedforward paths should be investigated as 


response-shaping methods. 
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APPENDIX A. COMPUTER PROGRAM INFORMATION 


The two computer programs are written in the common 
FORTRAN IV language for use on the IBM-360 model 67 computer 
With a CALCOMP'= model 765 plotter for graphics. This is the 
installation at the Naval Postgraduate School. Documenta- 
tion including data-card Puente information and internal 
secrvron Neadings is written into each program. Both programs 
are capable of expansion beyond the present limitations of 
tenth-order transfer function and third-order alpha/beta 
parameters with 16 curves per graph and 900 points per 


curve. Linear, semi-log and log-log plots are available. 


ite o> A let dn PROGR 

PARTF1 provides graphs of alpha, beta, and omega (ordi- 
nate) versus magnitude (abscissa). Numerator and denomi- 
nater polynomials in omega squared are formed. from the 
s-domain transfer function and magnitude squared determined 
from their quotient. Graph points are generated by incre- 
menting the desired ordinate parameter through the selected 
range of values while the other two parameters are held 
constant. The curve-family parameter is then set at the 
value for the next curve and the process repeated. Phase 
values are. determined along with magnitude whenever the 
incremented parameter is omega. These phase values are 


available for plotting along the abscissa against omega. 
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Subroutines used with PARTF1 are GRPLOT, ABVAR amd 
INTTFN. GRPLOT provides the plotting package for the 
program. The other two are used with the circuit simula- 
aber ie aanivire. 

ABVAR is used to enter the FORTRAN statements describ- 
ing the nonlinear parameter characteristics. Alpha and 
beta are initialized in the main program at each incremen- 
tation of omega, a solution made for magnitude and the 
information passed to this subroutine. The values of 
alpha, beta, magnitude and omega are then available for 
logic testing to determine the required alpha and beta 
values for that cirenit condition. Flag KNR will cetise a 
recycling of the magnitude solution using the new values 
of alpha and beta until the change in alpha is less than 
EPSA and in beta less than EPSB. KKl is available as a 
loop counter. It is initialized whenever alpha and beta 
ame initialized. 

In cases where one of the parameters is not a function 
Of the output transfer function, INTIFN is useay =e ae 
called by ABVAR when MDATA is set to one in the data set. 
The subroutine reads in additional data describing the 
internal transfer function and uses the ABVAR information 
to compute magnitude. This value, X2, is then passed back 
to ABVAR for use in determining parameter values. 

The storage requirements of PARTF1 are set primarily by 
the "DIMENSION PHASE (16,901), YY(16,901)" card. The number 


"16" represents the maximum number of curve parameters or 


i> 
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voltage input values used while "901" is the maximum number 
of computation points plus one. Using these values, the 
program requimees 205K DiS hor. Chew i.O0 tal mien 2.0 ame 

and about 19K for operations. This comes to about 233K. 
This may be reduced considerably by changing the statement 
mommetlcet The gach a lasanmbiege aint Hore ns Gancome 
mas six curves were to be plotted, the storage requirements 
would be reduced by about 72K (2 X 10 X 901 XK 4). If only 
500 points were used for each of the six curves, a further 
reduction in storage of about 19K would be gained. 

As expected, the number of curves and plotting points 
also affects the run time. Compile time is about 23 
seconds and the time for the link step about 1.5 seconds. 
ieee uke "ol thumb Wor. eo time is four seconds per 
900-point curve. A six-curve run of 900 points per curve 
would have an execute time of around 48.5 seconds or about 
50 seconds if a phase plot were also made. 


Figure 44 is a repeat of Figure 6 using only 300 points 


} 
per curve. Figure 45 uses only 100 points per curve. For 


design graphs of this size it would be wasteful in both 


storage and run time if 900 points per curve were demanded. 


2. PARTF2 PROGRAM 

‘-PARTFe provides graphs of alpha and beta versus omega, 
and alpha versus beta. As in the first program the input 
s-domain transfer function is separated into numerator and 
denominator functions which are individually transformed 


into the omega-squared domain. In PARTF2, however, the 
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denominator is then multiplied by the square of the magni 
tude and subtracted from the numerator to form avsadnglc 
polynomial in terms of the unknown ordinate parameter. The 
real roots of this polynomial together with the incremented 
abscissa values generate the graph points for the curve, 

The subroutines used with PARTF2 are GRPLOM. PORRT?> waa 
POLRT3. GRPLOT provides the plotting package for the 
program. The other two are root-finding subroutines. 
POLRT2 is used when the polynomial order is two and POLRT3 
when the order is three or more. Exponential overflow and 
underflow may occur in the high-order terms of the poly- 
nomial in POLRT3. This will not affect accuracy. “Te 
ERRSET command cancels the overflow and underflow crmonr 
print order and prevents the program from being prematurely 
ended when a given number of them have occurred. POLRT3 is 
a modified POLRT routine from the IBM System/360 Scientific 
Subroutine Package. 

Contrary to PARTF1 there is no one predominant factor 
setting the storage requirements. The limiting case of e@en 
factor contributes equally. Including the plotting package 
and operating space requirements, the storage requirement 
comes to about 148K. | 

PARTF2 run times are much longer than those for PARTFI1. 
Pormmng the polynomial, solving for the roovs, selecvuaee 
roots, etc. multiplies the run time over that for the simple 


quotient procedure of PARTF1 by a factor of seven. The go 


le 


~ 


time rule of thumb is 28 seconds per 900-point curve. The 
compile, link and edit time is about 23 seconds. Thus the 
estimated run time for eight 900-point curves would be 
abpouree4/7 seconds (6 X 28 4ez3). It pays to use fewer 


curves and fewer points with PARTFe. 
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